Entanglement Metrology Using a Joint Readout of Superconducting Qubits 
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Accurate and precise detection of multi-qubit entanglement is key for the experimental develop- 
ment of quantum computation. Traditionally, non-classical correlations between entangled qubits 
are measured by counting coincidences between single-shot readouts of individual qubits. We report 
entanglement metrology using a single detection channel with direct access to ensemble-averaged 
correlations between two superconducting qubits. Following validation and calibration of this joint 
readout, we demonstrate full quantum tomography on both separable and highly-entangled two- 
qubit states produced on demand. Using a subset of the measurements required for full tomography, 
we perform entanglement metrology with ~95% accuracy and ~98% precision despite ~10% fidelity 
of single measurements. For the highly entangled states, measured Clauser-Horne-Shimony-Holt 
operators reach a maximum value of 2.61 ± 0.04, and entanglement witnesses give a lower bound of 
~88% on concurrence. In its present form, this detector will be able to resolve future improvements 
in the production of two-qubit entanglement and is immediately extendable to 3 or 4 qubits. 



Since 1964, when Bell made testable [1] the famous 
paradox formulated by Einstein, Podolsky and Rosen [2] 
questioning non-locality in quantum mechanics, the mea- 
surement of correlations between quantum systems has 
been central to foundational tests against alternate theo- 
ries. The once 'spooky' non-classical correlations known 
as quantum entanglement have since been amply exam- 
ined in experiments [3, 4, 5, 6, 7]. Recently, entanglement 
has gained prominence as a key resource for large-scale 
quantum computation, making today's quantum engi- 
neer less concerned with foundational issues, and more 
with the generation and detection of near-perfect entan- 
glement between qubits. 

All approaches to metrology of entanglement, such as 
quantum state tomography, entanglement witnesses and 
generalized Bell inequality violations [8], require exper- 
imental measurement of ensemble-averaged qubit-qubit 
correlations. The paradigm for such detection used in 
trapped-ion systems [9, 10], where single-shot qubit read- 
outs with fidelities exceeding 99.99% have been realized 
[11], is to calculate correlations by measuring coinci- 
dences between individual detector 'clicks' over many 
repetitions. In solid-state systems, however, measure- 
ment cross-talk and lower readout fidelities make the ac- 
curate and precise measurement of correlations by this 
approach technically challenging, despite significant im- 
provements in recent years [12, 13, 14]. Various ap- 
proaches seeking direct access to correlations with joint, 
or quadratic, on-chip detectors have been theoretically 
proposed for mesoscopic systems [15, 16]. 

In this Report, we demonstrate entanglement metrol- 
ogy with a single measurement channel that performs 
a joint readout of two charge-based superconducting 
qubits. This readout is based on qubit-state depen- 
dent shifts of a microwave transmission-line resonator 



in a circuit quantum electrodynamics architecture [17]. 
By validating and calibrating the measurement model 
for this single-channel detector, we circumvent a low 
~10% single-shot detection fidelity and obtain ensemble- 
averaged correlations with ~95% accuracy and ~98% 
precision. While this joint readout has been previously 
used in superconducting qubit systems [18], we demon- 
strate for the first time its use to detect state-of-the- 
art highly-entangled states produced on demand, with 
fidelity to targeted Bell states >90%, a lower bound of 
~88% on concurrence, and violation of Clauscr-Horne- 
Shimony-Holt (CHSH) inequalities [19] with a value of 
2.61 ± 0.04. Furthermore, we anticipate that future im- 
provements in entanglement production with this archi- 
tecture will be detectable with the joint readout in its 
present form. 

Our device consists of a superconducting transmission- 
line resonator that couples two transmon qubits [20, 21]. 
Beyond mediating the interaction between the qubits, 
the bus serves as a single detection channel for their 
joint readout. Arbitrary single-qubit x- and y-rotations 
are performed using in-phase and quadrature microwave 
pulses resonant with the ground to first-excited state 
transition of each transmon [22]. Reduction of leakage 
to the second-excited state is accomplished through the 
technique of derivative removal by adiabatic gate [23], 
resulting in error rates of ~1% in single-qubit rotations 
(See attached supporting material). Local flux-bias lines 
tune the qubit transition frequencies on nanosecond time 
scales, allowing control of single-qubit dynamical phases 
and of a Z h t$Z n interaction, both crucial for implement- 
ing a two-qubit conditional-phase (C-phase) gate [24] 
(here Z h ^ = cr^ R ^ is the single-qubit Pauli z opera- 
tor [25] on the left (right) qubit). A schematic of the 
experimental setup [26] is shown in Fig. 1A. 
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FIG. 1: Experimental setup and readout. (A) Circuit diagram of experimental setup, representing the cavity as a lumped L — C 
resonator. Microwave drive tones that address the cavity, lorf, and qubits, <^' R \ are applied via the cavity input line. The left 
(right) qubit has charging energy E^,^ and Josephson energy _E^ R ' (flux $ L ( R ' tunable via y L ( R >). The multiplexed qubit- 
state information is transmitted out on a single readout line, amplified through a high electron mobility transistor (HEMT) 
amplifier at 4K and mixed down at room temperature for digital processing. (B-E) Color images of transmitted Vh amplitude 
as a function of time and cavity drive frequency wrp, for two-qubit states (B) |0, 0), (C) jO, I), (D) |1, 0), and (E) |1, 1). In 
this color scheme blue (white) represents Vh = 0(1.9)mV. On short time scales t < Tj , maximal transmission occurs at 
different frequencies for the four cases (see text for details). (F) Temporal average of the homodyne voltage transients in (B-E) 
over the first 500 ns shows well resolved peaks. (G) Measured Vh transients for drive at the cavity frequency corresponding to 
|0, 0), showing the measurement transient for the four computational states. 



Combining the single-qubit rotations with the C-phase 
gate, it is possible to produce maximally-entangled states 
such as the four standard Bell states, \$ ±) = (|0,0) ± 
|l,l))A/2 and |$±) = (|1,0) ± |0,l))A/2, where |/,r) 
denotes excitation level l(r) of the left (right) qubit [24]. 
Yet, to accurately and precisely detect arbitrary two- 
qubit states, a complete physical model and calibration 
of the joint readout is necessary. 

The physical mechanism enabling the joint readout is a 
qubit state-dependent dispersive cavity shift that is large 
relative to the cavity linewidth k. In this 'strong disper- 
sive' regime [27], the system is described by a dispersive 
Tavis-Cummings Hamiltonian 

H TC /h = ( Wc + X L Z L + X R Z R ) a t a - - ^ R ,(l) 

where wc is the bare resonator frequency, w L ( R ) is the 
Lamb-shifted ground to first excited state transition fre- 
quency for the left (right) qubit, and x L is the left 



(right) qubit-state dependent cavity shift. Actual param- 
eter values were determined by spectroscopy experiments 
[26]. From Eq. (1), there can be at most four distinct cav- 
ity frequencies corresponding to the joint state of the two 
qubits. 

The linear dispersive shifts are calibrated by mea- 
surement of the transmitted homodyne voltage Vh 
transient when pulsing a measurement tone, hav- 
ing prepared each of the four computational basis 
states, |0, 0), |0, 1), |1, 0), |1, 1), with single-qubit rota- 
tions. Color images in Figs. IB-IE show the transient 
(Vh) as a function of cavity drive frequency wrf for the 
four cases (brackets denote an average over 600,000 rep- 
etitions). On time scales t < T^ K ' = 1.2(0.9) /js, the 
largest transmission occurs at distinct frequencies for all 
four cases. The discrete transmission peaks are well re- 
solved in Fig. IF, showing the time average of (Vh) versus 
frequency over the first 500 ns. By matching these fre- 
quencies to Eq. (1), we extract x L(R) / 27r = 13(4) MHz. 
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On long time scales t > , the dominant transmis- 

sion is at the peak corresponding to |0, 0) for all cases. 
The transition between the two time limits is most evi- 
dent for 1 1 , 1 ) , which decays partially into |0, 1) and |1, 0) 
before completely relaxing into the ground state |0, 0). 

It is because % L and x R are both larger than k = 
1 MHz that a measurement can give joint qubit informa- 
tion. For example, applying a measurement tone at the 
cavity transmission peak corresponding to 1 , 0) queries 
for the joint property that both qubits are in their ground 
state. Given the large state-dependent dispersive shifts, 
the system being in any other state results in a low trans- 
mission signal. This is best shown by the four transmis- 
sion transients in Fig. 1G. The transients for the states 
|0, 1), |1,0) and |1,1) are all different from the transient 
for |0, 0), but also not identical to each other. Qubit re- 
laxation during the measurement results in the transients 
converging towards the |0, 0) response on long time scales. 

More rigorously, the idealized measurement M is a pro- 
jection operator onto |0,0), M = |0,0) (0,0| = (I + Z L + 
Z R + Z L &Z R )/4, which is sensitive to the polarization of 
each qubit along its z-axis as well as to two-qubit corre- 
lations. However, the actual measurement is an ensemble 
average of Vh which, due to qubit relaxation and partial 
overlap of the dispersive peaks, we expect to be described 

by 

(Vh) = Pn + Pzi(Z h ) + piz{Z K ) + Pzz{Z h ® Z K ) + 5v, 

(2) 

where the 0lr are constant coefficients and Sv is classical 
amplifier noise (details in supporting material) . The am- 
plifier noise, which limits the single-shot fidelity to ~10%, 
is largely eliminated by repeating the state preparation 
and measurement 600,000 times over fifteen seconds. 

A simple set of Rabi flopping experiments can be per- 
formed (Figs. 2A-C) to validate the measurement model 
of Eq. (2) and place bounds on other terms. The most 
general model, including all linear combinations of two- 
qubit Pauli operators [28], is 

M= PlrL®R. (3) 

L,Re{I,X,Y,Z} 

Figures 2A and 2B show the measured (Vh) as a func- 
tion of the duration of an applied drive at ui L and ui R , 
respectively. In each case the drive induces a Rabi os- 
cillation of the addressed qubit around the y-axis of its 
Bloch sphere [25]. The observed oscillations in (Vh) are 
fit (solid lines) excellently by assuming the measurement 
model in Eq. (3) and the theoretical (Z) and (X) for 
driven qubits (see details in supporting material). From 
the fit to Fig. 2A (2B), we can estimate deviations from 
Eq. (2) due to the terms (X L ) ((X R )) and (XZ) ((ZX)) 
are each < 2% of the full range of (Vh)- 

The third experiment (Fig. 2C) measures the homo- 
dyne response to simultaneous Rabi drives on both 
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FIG. 2: Measurement model from Rabi oscillations. (A-C) 
Rabi oscillations on the (A) left qubit, (B) right qubit, and 
(C) simultaneously on both. Solid lines are fits to the model 
in Eq. (3). See text for details. (D) Fourier transform (FT) of 
the three Rabi experiments (markers) and of best fits (curves). 
While the red (blue) traces show one main peak at the Rabi 
frequency f2 L( - R ', the purple traces reveal peaks at fi L , f2 R , 
fi L + 17 R , and f2 L — f2 R , demonstrating the mixing property 
that makes the joint measurement sensitive to qubit-qubit 
correlations. 



qubits. The observed oscillations not only show fre- 
quency components at the individual Rabi frequencies, 
f2 L and f2 R , but also at their sum and difference. This 
is clearly revealed in the Fourier transform of the oscilla- 
tions (Fig. 2D). This mixing effect agrees quantitatively 
with the term (ZZ) in Eq. (2). An excellent fit is also 
obtained, giving deviations from Eq. (2) due to {XX) 
that are ~ 2% of the full swing of (Vh) (see details in 
the supporting material). These Rabi experiments thus 
corroborate the measurement model Eq. (2) . 

Besides testing the measurement model, these Rabi- 
flopping experiments allow calibration of the coefficients 
P LR in Eq. (2). The best fits give (3 n = 800 fiV, f3 IZ = 
380 /iV, (3 Z i = 380 ^V, and (3 Z z = 200 fiV. The large 
ratios Pzz / Piz(zi) ~ 0-6 indicate the high sensitivity of 
the joint readout to the qubit-qubit correlations. 

Having substantiated the physical nature and the 
quantitative model of the joint readout, we next use it to 
perform tomography of a variety of two-qubit states, both 
separable and entangled, produced on demand. We first 
perform an over-complete set of 30 raw measurements, 
each obtained by applying a different pair of simultaneous 
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FIG. 3: Representation of two-qubit states using the Pauli set. 
Experimental Pauli set (with trivial (II) = 1 not shown), ob- 
tained from linear operations on raw measurement data, for 
separable states (A) |1,0> and (B) (|0,0) - |0, 1) + |1,0) - 
1, l))/2 and entangled states (C) |*+) and (D) the Bell state 
in the z-basis \$+) x = (|0, 0) + |0, 1) - |1, 0) + |1, l))/2. Red 
(blue) bars correspond to left (right) single-qubit Pauli opera- 
tors. Purple bars are the qubit-qubit correlations. The fideli- 
ties to the four ideal targeted states computed using Eq. (5) 
are T = 98.2±0.4%, 96.8±0.4%, 90.0±0.6%, and 92.5±0.7%. 
The density matrix obtained using the same raw measure- 
ments is shown in three-dimensional city-scape format to the 
right of each Pauli set. 



rotations of the qubits prior to measurement. The qubit 
rotations are chosen from the set {/, R±\R^ /2 ,R^ /2 } 
(see supporting material). Using the calibration of the 
Plr and least-squares linear estimation, we then con- 
struct the Pauli set P, whose 16 elements are the ex- 
pectation values of the two-qubit Pauli operators, (LR), 
where L,R g {I, X,Y, Z}. Two-qubit states can then 
be visualized by plotting the Pauli set in bar-graph for- 
mat, in contrast to plotting the density-matrix in three- 
dimensional city-scape format [25]. 

An advantage of the Pauli set representation is that 
one can easily distinguish separable from entangled 
states. In Fig. 3, we use the Pauli set to view the ex- 
perimentally generated separable states, (A) |0, 1) and 
(B) (|0) + |1)) L ® (|0) - |1))r/2, as well as entangled 
states, (C) |4>+), and (D) the Bell state in the x-basis 
|*+) a = (|0,0) + |1,0) - |0,1> + |l,l»/2. For these 
four states, the Pauli set ideally contains three non- 



zero bars, all of unit magnitude. The Pauli set can 
be sub-divided into three sections: the qubit polar- 
ization vectors, P L = {(X L ), (Y L ), (Z h )} and P R = 
{(X R ), (Y R ), (Z R )}, and the vector of two-qubit corre- 
lations Q = {(XR),- ■ ■ , (ZZ)}. In the figure, P L , P R 
and Q are color-coded red, blue and purple, respectively. 
For the separable states, we observe near unity compo- 
nents in the three sub-sections of the Pauli set, P L , P R , 
and Q. In contrast, for the entangled states, the only 
large components are in Q. The presence of large bars in 
Q and vanishing P L , P R is a direct signature of a high 
degree of entanglement. 

The Pauli set representation permits testing of some 
simple physical constraints, which if not met, reveal sys- 
tematic errors. The most easily tested physical constraint 
is \(LR)\ < 1, \/L,R (other constraints are < |P L ' R | < 
1, and < |Q| < Vo). Although the states shown in 
Fig. 3 satisfy these bounds, a more thorough experimen- 
tal test for systematic errors is to measure the Pauli set 
for a collection of states that differ only by the angle of 
a single-qubit rotation prior to measurement. Two such 
evolutions are shown in Figs. 4A and 4B, which involve a 
rotation 6 of the left qubit about its y-axis after preparing 
the separable state |0, 0) (evolution A) and the entangled 
state (evolution B), respectively. 

In evolution A, varying 6 produces the oscillation of the 
Pauli set shown in Fig. 4C. Systematic errors in detection 
could appear as offsets and amplitudes of the Pauli set 
betraying the ±1 bounds. Such deviations would limit 
the accuracy of physical quantities extracted from the 
Pauli set, and thus are important to identify and correct. 
In Fig.4C, (XI), (Zl), (XZ), and (ZZ) oscillate with 
an average visibility of 97.6 ± 0.3%, demonstrating the 
large swing of the meter. Moreover, the average absolute 
error of all the ideally zero- valued bars is < 10%. In evo- 
lution B, the dominant oscillating components are all in 
Q, indicating that the state remains entangled through- 
out all the rotations. In this case, we find a visibility 
of 91.5 ± 0.3%, in good agreement with a master equa- 
tion simulation taking into account qubit relaxation and 
dephasing. A residual oscillation amplitude of ~10% is 
observed in (XI) and (Zl), a factor ~2 larger than ex- 
pected from theory. This discrepancy can arise from a 
combination of small calibration errors in single-qubit ro- 
tations and various residual higher order couplings (see 
supporting material). In the future, composite pulses 
and shaping might further reduce these effects. 

Finally, we use the fully characterized and high- 
visibility readout to extract accurate and precise mea- 
surements of physical quantities intrinsic to the two-qubit 
state. Examples include the fidelity to the targeted state, 
entanglement witnesses, and CHSH operators. Interroga- 
tions of such quantities are henceforth statements about 
the quality of the states prepared. 

The fidelity T to a targeted state provides one way of 
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FIG. 4: Evolution of the Pauli set and entanglement witnesses under qubit rotations. (A and B) Gate sequences generating 
evolutions of the Pauli set by subjecting (A) the separable state |0, 0} and (B) the entangled state to a rotation R y (9) on 

the left qubit, —210° < 9 < 210°. (C and D) The measured Pauli set as a function of 9 for evolutions A and B, respectively. 
In evolution A, the left qubit polarization rotates along the x-z plane, while the right qubit remains fully-polarized along z. In 
evolution B, both qubit polarizations vanish, with the only nonzero and oscillating Pauli operators being qubit-qubit correlators 
(purple bars). Arrows at 9 = — 90(+90)° indicate when the ideal two-qubit state is the standard Bell state |3>-) (|^+)). (E 
and F) Experimental lower bounds Bi (orange) on the concurrence given by the optimal witnesses for Bell states W*+ (circles), 
W*_ (squares), W*+ (triangles), and W#- (crosses), and fidelity T to the ideal state (black circles) for (E) evolution A and (F) 
evolution B. In F, a maximum lower bound is reached by B<s>- (£?*+) at 9 — — 90(+90)°. (G-H) Experimental average value of 
CHSH operators Czzxx (circles), Czxxz (squares), Czxzx (triangles), Cxxzz (crosses). (G) For evolution A all (C) values 
stay within the separable state bounds ±v2 up to measurement noise, while for (H) evolution B, max|(C)| = 2.61 ± 0.04. 
Solid lines are master equation simulations. 



quantifying the control over two-qubit states and is given 
by the inner (dot) product of the measured P to that of 
the ideal state, 

T = -P ■ Pidoal- (4) 

For the states produced in evolutions A and B, we find 
average fidelities T = 98.8±1.0% and 93.4±1.5% over all 
8, respectively. The fidelities (black circles) in Figs. 4E-F 
show excellent agreement with simulation (solid lines). 
This agreement demonstrates the accuracy of our meter 



for detecting both separable and entangled states, and 
the high fidelity values attest to the level of two-qubit 
control in the experiment . 

Measures beyond fidelity are necessary to quantify the 
degree of two-qubit entanglement. Often, entanglement 
monotones such as concurrence C [8] are obtained us- 
ing non-linear estimators. It is standard to first perform 
maximum-likelihood estimation [29] to generate a phys- 
ical p despite any statistical or systematic errors in the 
raw data, and to then calculate these metrics from the 
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eigenvalue spectrum of related matrices [8]. This non- 
linear process complicates the propagation of any statis- 
tical and systematic errors in the measurements, and can 
bias the estimation of such metrics as the purity of the 
two-qubit state increases [30, 31] (See supporting mate- 
rial). 

To be quantitative about entanglement while using 
only linear operations on the raw measurements, we 
make use of entanglement witnesses [8, 32]. An entan- 
glement witness is a unity-trace observable W with a 
positive expectation value for all separable states, such 
that Tr(pW) < guarantees entanglement. Further- 
more, — 2Tr(pW) gives a lower bound on C [32]. The 
optimal witness for a targeted entangled state gives the 
tightest lower bound. For the Bell states, these are 

W*+ = ^(11 - XX + YY - ZZ), 

Wv- = 1 -{II + XX-YY -ZZ), 

W$+ = - XX - YY + ZZ), 

W$_ = 1 -{II + XX + YY + ZZ), 

giving lower bounds Bi — — 2Tr(pWi) on C. 

The four bounds obtained for evolution A(B) are 
shown in Fig. 4E(4F). In evolution A, the four bounds 
arc non-positive for all 8 to within measurement error, 
indicating that entanglement is not witnessed. This is 
as expected, since single-qubit rotations should not pro- 
duce any entanglement. In Fig.4F, in contrast, bounds 
B<s, + and extend into the positive region, reaching 
85.9 ±1.5% and 88.1 ±1.5% at 8 = -90° and 90°, respec- 
tively. There is at least one positive bound for most 8 
(excluding ±180° and 0°), indicating that the two qubits 
are entangled. Solid lines in Figs. 4E-F are obtained from 
the master equation simulation. The agreement between 
data and theory shows the accuracy of the entanglement 
witnesses and the small residuals ~2% demonstrate their 
precision. 

A more well-known entanglement measure is the CHSH 
operator, often used to test local-hidden variable theo- 
ries. CHSH operators [19] are most generally defined as 

C AtB ,A',B'=AB + AB' + A'B-A'B', (5) 

with left-qubit operators A, A' and right-qubit operators 
B, B' being single-qubit Pauli operators along any axis. 
For separable states, |(C)| < 2. However, for A J- A' and 
B _L B', the separable bound is tighter, [(C) | < \f2. 

From the Pauli sets of evolutions A and B, we ob- 
tain experimental expectation values of four CHSH op- 
erators with A, A' e {X L ,Z L } and B,B' e {X R ,Z R }. 
For the separable states of evolution A (Fig. 4G), we find 
that to within statistical error, (Czxxz) = {XX) — 
(XZ) + (ZX) + (ZZ) (squares) and (Czxzx) = (XX) + 



(XZ) — (ZX) + (ZZ) (triangles) remain within the sep- 
arable bounds for all 9. For the entangled states pre- 
pared in evolution B (Fig. 4H), instead, (C) oscillates 
past the separable bounds. At 8 — ±45°, a maximum 
value |(C) | = 2.61 ±0.04 is reached. The agreement with 
theory and proximity of the maximum |(C) | to the 2y/2 
upper bound [33] further demonstrate the high degree of 
entanglement of the states produced. We note that while 
(C) exceeds the local-hidden variable bound [8] of 2 by 
~15 standard deviations, this result is not a strict viola- 
tion of local-hidden variable theories due to locality and 
measurement loopholes in our system. For superconduct- 
ing qubits, the measurement loophole has recently been 
closed using high-fidelity single-shot readout [7]. 

In summary, we have performed metrology of en- 
tanglement between two superconducting qubits using 
a single measurement channel giving direct access to 
qubit correlations. While loopholes preclude using 
this joint readout for fundamental tests of quantum 
mechanics, the demonstrated ability to measure fidelity 
to targeted states, lower bounds on concurrence and 
strong violation of CHSH inequalities is eminently useful 
for quantum-computer engineering. In its present form, 
this single-channel detector will be able to detect future 
improvements in the controlled production of entangled 
states and is immediately extendable to three or four 
qubits. Similar implementations of joint detection may 
prove useful for metrology of multi-qubit operations 
in other qubit platforms where cross-talk and low 
measurement fidelity make single-shot individual qubit 
readouts challenging. 
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I. MATERIALS AND METHODS 
A. Sample Fabrication 

The device is fabricated on a 430 ^im thick sapphire 
substrate. The superconducting transmission-line cav- 
ity and flux-bias lines are defined via optical lithography 
and fluorine-based reactive ion etching of a dc-sputtered 
niobium film (180 nm thick). The two transmons are 
patterned using electron-beam lithography with split- 
junctions, grown using double-angle deposition of alu- 
minum, with layer thicknesses of 20 nm and 90 nm. The 
sample is cooled in a dilution refrigerator to 13 mK. The 
experimental setup is presented in detail in a previous 
publication [1]. 

B. Experimentally Determined Parameters 

From a set of heterodyne transmission measurements 
obtained when tuning each qubit into near resonance 
with the cavity, we determine the qubit-cavity coupling 
strengths 5 L (r)/2tt = 199(183) MHz. Fitting a multi- 
level Tavis-Cummings Hamitonian to spectroscopy mea- 
surements of the two lowest transitions of each qubit, 
we extract maximum Joscphson energies E™^ K -./h = 
28.48 (42.34) GHz and electrostatic charging energies 
E CMK)/h = 317(297) MHz. Standard sliding 7r-pulsc 
and Ramsey fringe experiments give qubit relaxation 
times T^ R ) = 1.2 (0.9) and dephasing times — 
1.5 (1.1) /js. 



II. EXPERIMENTAL REDUCTION OF HIGHER 
LEVEL LEAKAGE 

For superconducting qubits, previous work [2, 3] has 
demonstrated single-qubit gate errors of 1-2%. Specifi- 
cally, in Ref. [3], this average gate error is consistent with 
qubit relaxation and dephasing for gate lengths t ga tc > 16 
ns, and with leakage to the second excited level at shorter 
gate lengths. Recent theory proposes a possible reduc- 
tion of errors due to this leakage to the second excited 
level via the technique of derivative removal via adiabatic 



gate (DRAG) [4]. In DRAG, when a single-qubit rota- 
tion pulse is applied along the derivative of the 
pulse is applied along the y-axis to cancel out the higher 
level leakage, and vice versa. 

In this work, we have used DRAG on all single-qubit 
gates. Standard x- and y-rotations are performed with 
in-phase and quadrature microwaves tuned to the qubit 
ground to first excited state frequency, and are shaped 
with Gaussian envelopes, truncated to two standard devi- 
ations a on each side. After each gate, a 5 ns buffer is in- 
cluded to avoid any overlap with the following gate. The 
leakage due to an x (y)-rotation or in-phase (quadrature) 
pulse is reduced by applying a complementary quadra- 
ture (in-phase) tone shaped by a truncated derivative-of- 
Gaussian envelope. 

Using the technique of randomized benchmarking (RB) 
[5], we extend the extraction of gate errors as performed 
in [3] to our single-qubit gates with and without DRAG, 
finding a significant improvement in the error rates with 
DRAG. As shown in Fig. SI, we find that without DRAG 
(red circles), the gate errors decrease with decreasing 
gate length before turning back up at a gate length of 
17 ns, signaling the onset of second-excitation leakage. 
The overall higher and earlier onset of the increase in 
gate error compared with Ref. [3] is due in part to the 
shorter T\ and the larger g of the present qubit. In 
contrast, with DRAG (blue circles), the gate errors con- 
tinue to decrease, reaching 0.7±0.5% at the shortest gate 
lengths possible with our arbitrary waveform generator 
(Tektronix 5014B). 

The total gate length chosen for this work is 13 ns 
((7 = 2 ns) at which the error per gate is 0.9 ± 0.5%, in 
excellent agreement with the simple theory (solid black 
line) with T\ = 1.2 ^s for the left qubit. In the future, 
using DRAG with improved pulse-shaping resolution in 
hardware and improved coherence times could further 
lower single-qubit error rates towards the fault tolerant 
threshold [6]. 

III. MEASUREMENT MODEL VALIDATION 

As mentioned in the main text, we use fits to the three 
Rabi-flopping experiments of Fig. 2 to place bounds on 
deviations from the measurement model of Eq. (2). Be- 
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evolution of (Z) and (X) is given by 
7i72 + i2^ 7i72 + " 




si- 71721-"- 

- [2 1 l 2 + 7l (7 2 -7 1 )]sin(fii)^ 

cos(S2t) — = 

20 



Here, Si = ,/Sl 2 - (1/t r ) 2 is an effective oscillation rate, 
7i = l/Ti is the relaxation rate, 72 = 71 /2 + 7^ is the 
dephasing rate, and tr — 2/(71 + 72) is the Rabi decay 
time. 

Using these expressions in Eq. 3 and fitting to the three 
experiments, we can estimate the coefficients 0lr. For 
single-qubit driving [Fig. 2A(B)], the right (left) qubit is 
always in the ground state, and only terms (Z L ), (X L ), 
(XZ) and (ZZ) ((Z R ), (A R ), (ZX), and (ZZ)) con- 
tribute to the (Vh) oscillation. Fitting the form 

(Vb)a<b) =W + W x {Z^) + W 2 (X L W), 



FIG. SI: Reduction of single-qubit gate errors via 

DRAG. Error per gate for the left qubit extracted from ran- 
domized benchmarking for different gate lengths using both 
standard Gaussian-shaped pulses (red squares) and DRAG- 
enhanced pulses (blue squares). Excellent overlap with a 
model for gate error including qubit relaxation (black curve) 
suggests that DRAG-enhanced pulses successfully eliminate 
leakage to the second-excited state. Gate errors ~ 1%, which 
are otherwise unattainable with standard pulses, are reached 
using DRAG. 



L(R) 

7i \ 



with W , W u W 2 , S1 L(R) , 
rameters gives an excellent fit. In both cases, the best-fit 



j L(R) r 

and 7 2 as free pa- 



is less than 



W 2 , corresponding to X i{ix) + 0xz(zx), 
2% of the full range of (Vh), ~ Wiz + 20zi- For the 
doubly-driven case (Fig. 2C), the fit function used is 

(V H )c= 0ii + Pxi(X l ) + 0zi (Z L )+ 0i x (X R ) + IZ (Z R ) 
+0xx(XX) + 0xz(XZ) + 0zx(ZX) + 0zz(ZZ), 

with 0ij, S1 L , S1 R , 7 L , and 7J 1 as fit parameters. The 
best-fit coefficients captured in Eq. (2) are 0n — 800 fiV, 
0iz = 380 /iV, 0zi = 380 fiV, Z z = 200 ^V. Best-fit 
values of the remaining coefficients are each less than 2% 
of the full range of (Vh)- 



cause in these tests each qubit is driven around the y- 
axis of its Bloch sphere, all terms involving Y h and Y R 
in Eq. 3 would not contribute to (Vh)- The presence of 
such terms can be tested by rotating each or both qubits 
around their x-axis instead. We do not find any signif- 
icant differences in such experiments from the ones pre- 
sented in the text, and the results here can be generalized 
for both quadratures X and Y. 



In our experiment the detuning ^1.5 GHz between the 
two qubits is large compared to the Rabi-flopping rates, 
and we can assume a simple model of independent qubit 
driving. For a qubit driven at a rate SI around its y- 
axis starting from the ground sate, the theoretical time 



IV. TWO-QUBIT STATE TOMOGRAPHY 

Full tomography of the two-qubit state is performed 
by using an over-complete set of 30 raw measurements. 
These measurements involve applying different simulta- 
neous rotations on the qubits, as listed in Table SI. 
The 15 measurements labeled P t involve positive rota- 
tions chosen from {I, R+" , Fit*' 2 , Rf /2 }. The remain- 
ing 15, labeled iVj, involve negative rotations chosen from 
{I,R- n ,R^ 7T/2 ,Ry' K/2 }. Ensemble averages of Pi and iVj 
are obtained by repeating state preparation, analysis ro- 
tation, and measurement over 600,000 times. A standard 
least-squares linear estimator is then used to extract the 
Pauli set P discussed in the text. Although just 15 lin- 
early independent measurements (such as either all Pi 



3 



or all Ni) is necessary for state tomography, using all of 
these rotations and least-squares estimation reduces the 
statistical error in the extraction of P. 
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TABLE SI: The 30 raw measurements. 



FIG. S2: Biasing of entanglement bounds by MLE. 

Comparison of the lower bound £?,/,_ on concurrence C com- 
puted from simulated noisy raw data with and without use 
of maximum-likelihood estimation. The bound computed us- 
ing MLE systematically underestimates the true bound (in 
this case, always equal to the true concurrence C, red line), 
while the bound computed directly from the simulated raw 
data remains faithful even as the Werner state approaches 
the pure Bell state i.e., as A — > 1. For a pure Bell state, 

the MLE-computed bound underestimates the true bound by 
4%. 



sets of simulated raw measurements for each A by assum- 
ing Gaussian amplifier noise consistent with the experi- 
ment. Figure S2 shows the concurrence bound £>^_ as a 
function of the true C of the Werner state, obtained with 
and without MLE processing of the simulated noisy data. 
We find that while the mean of estimated directly 
from the raw data is unbiased, the mean of the concur- 
rence bound obtained with MLE becomes increasingly 
biased the more pure the Werner state, i.e., the closer 
A is to unity. MLE underestimates the bound by 1% at 
C = 0.85, and by 4% at C = 1. 



V. BIASING OF METRICS BY 
MAXIMUM-LIKELIHOOD ESTIMATION 

Maximum- likelihood estimators (MLE) can become bi- 
ased if the true mean lies close to a boundary of the 
allowed parameter space. In order to quantify the im- 
portance of this effect on the estimation of lower bounds 
on concurrence C given by entanglement witnesses, we 
have performed Monte-Carlo simulations for nearly-pure 
Werner states [7], 

p w (X) = \\ip-)(ip-\ + (l-\)I/4, 

with Werner parameter A € [0.8, 1]. We have created 100 



VI. POSSIBLE SYSTEMATIC ERRORS 

There are a variety of higher order and systematic ef- 
fects that affect the accuracy of the entanglement metrol- 
ogy at the few-percent level. For example, the discrepan- 
cies between the experiment and master equation simula- 
tion in Fig. 4E can arise from a systematic under-rotation 
of both qubits by only 1%. There are also higher order 
couplings that are relevant at this level. The first is the 
finite strength of the two-qubit ZZ entangling interaction 
even in the off state (C/2tt ~ 1.2 MHz [1]). This residual 
coupling leads to errors in some of the two-qubit correla- 
tors on the order of £/f2 L ( R ) ~ 2%. A second is the pres- 
ence of the qubit-qubit swap interaction (J/2ir ~ 60 MHz 



[8]), leading to errors of order J/(w L — w R ) ~ 4%. An- 
other effect is the qubit-state dependent filtering of the 
drive applied to a qubit, which is expected to be on the 
order of X R ( L )/( W L ( R ) - Wc ) ~ 2%. The effect of these 
couplings can be mitigated by implementing appropriate 
composite pulse schemes [9] and will be explored in the 
future. 



[1] L. DiCarlo, et al, Nature 460, 240 (2009). 

[2] E. Lucero, et al, Phys. Rev. Lett. 100, 247001 (2008). 

[3] J. M. Chow, et al, Phys. Rev. Lett. 102, 090502 (2009). 

[4] F. Motzoi, J. M. Gambetta, P. Rebentrost, F. K. Wilhelm, 
arXiv.0901.0534 (2009). 

[5] E. Knill, et al, Phys. Rev. A 77, 012307 (2008). 

[6] D. Gottesman, Stabilizer codes and quantum error cor- 
rection, Ph.D. thesis, California Institute of Technology, 
Pasadena (1997). 

[7] R. Horodecki, P. Horodecki, M. Horodecki, K. Horodecki, 
Rev. Mod. Phys. 81, 865 (2009). 

[8] J. Major, et al, Nature 449, 443 (2007). 

[9] H. K. Cummins, G. Llewellyn, J. A. Jones, Phys. Rev. A 
67, 042308 (2003). 



